Next-to-leading BFKL phenomenology of forward-jet cross sections at HERA 
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We show that the forward-jet measurements performed at HERA allow for a detailed study of 
corrections due to next-to-leading logarithms (NLL) in the Balitsky-Fadin-Kuraev-Lipatov (BFKL) 
approach. While the description of the da/dx data shows small sensitivity to NLL-BFKL corrections, 
these can be tested by the triple differential cross section da / dxdk\dQ 2 recently measured. These 
data can be successfully described using a renormalization-group improved NLL kernel while the 
standard next-to-leading-order QCD or leading-logarithm BFKL approaches fail to describe the 
same data in the whole kinematic range. We present a detailed analysis of the NLL scheme and 
' renormalization-scale dependences and also discuss the photon impact factors. 
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^\ I. INTRODUCTION 

Forward-jet production in lepton-proton deep inelastic scattering is a process in which a jet is detected at forward 
rapidities in the direction of the proton. This process is characterized by two hard scales: Q 2 , the virtuality of the 
\ intermediate photon that undergoes the hadronic interaction and fcf,, the squared transverse momentum of the forward 
jet. When the total energy of the photon-proton collision W is sufficiently large, corresponding to a small value of 
y—( \ the Bjorken variable x~Q 2 /W 2 , forward-jet production is relevant [l| for testing the Balitsky-Fadin-Kuraev-Lipatov 
\Q ■ (BFKL) approach 0. 

In fixed-order perturbative QCD calculations, the hard cross section is computed at fixed order with respect to 
a s , and large logarithms coming from the strong ordering between the proton scale and the forward-jet scale are 
resummed using the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equation However in the 
small— a; regime, other large logarithms arise in the hard cross section itself, due to the strong ordering between the 
energy W and the hard scales. These can be resummed using the BFKL equation, at leading (LL) and next-leading 
i-^ ' (NLL) logarithmic accuracy 0, H[ . 

It has been shown that the HI and ZEUS da/dx forward-jet data [1,[1,0| are well described by LL-BFKL predictions 
Q_i! [1, El, while fixed-order perturbative QCD predictions at next-to-leading order (NLOQCD) fail to describe the data, 
underestimating the cross section by a factor of about 2 at small values of x. However, these tests on the relevance 
of BFKL dynamics have not been considered fully conclusive. On the theoretical side, it has been found that NLL- 
BFKL corrections Q could be large enough to invalidate the tests. On a phenomenological side, other models such 
as DGLAP evolution with a "resolved" photon [l(| could increase the NLOQCD predictions and come to reasonable 
agreement with the data. 

The recent experimental forward-jet measurements [H, Q performed at HERA motivate a new phenomenological 
analysis of BFKL effects in forward-jet cross sections. In particular the triple differential cross section da / dxdk^dQ 2 , 
allows for a detailed study of the QCD dynamics of forward jets. Contrary to the da/dx data, which were obtained 
with kinematical cuts such that r — k^/Q 2 ^ 1, the triple differential cross section is measured with different sets of 
cuts such that the data are also sensitive to the regime r ^ 1 , where the two hard scales of the problem are somewhat 
ordered. While LL-BFKL predictions describe well the data obtained with r~ 1, it was noticed [§] that they fail to 
describe the r3>l regime, indicating the need for NLL-BFKL corrections. 
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It was known that NLL-BFKL corrections could be large due to the appearance of sp urious singularities in contra- 
diction with renormalization-group requirements. However it has been realized (TH Il2j | that a renormalization-group 
improved NLL-BFKL regularization can solve the singularity problem and lead to reasonable NLL-BFKL kernels 
(see also [l3[ for different approaches). This motivates the present phenomenological study of NLL-BFKL effects in 
forward-jet production. Even though the determination of the next-leading impact factors is still in progress [ijj], 
our analysis allows us to study the NLL-BFKL framework, and the remaining ambiguity corresponding to the de- 
pendence on the specific regularization scheme. Our goal is to confront the new experimental data, in particular the 
triple-differential cross section, to NLL-BFKL predictions in different schemes. 

In Ref . [l5j . such a phenomenological investigation has been devoted to the proton structure function data, taking 
into account NLL-BFKL effects through an "effective kernel" (introduced in pjj) using different schemes. A saddle- 
point approximation for hard enough scales was used to evaluate the BFKL Mellin integration which allowed one to 
obtain a phenomenological description of NLL-BFKL effects. In the present study devoted to forward-jet production, 
we take into account the proper symmetric two-hard-scale feature of the forward-jet problem when introducing the 
effective kernel, and we implement the NLL-BFKL effects with an exact Mellin integration, rather than a saddle-point 
approximation. 

Some preliminary results, mostly based on the saddle-point approach, were presented in (l6| . They showed the 
potential of forward-jet data on da/dx and specially da / dxdk^dQ 2 to discuss NLL effects in the BFKL approach. In 
this paper, we systematically use an exact Mellin integration, and we present a detailed analysis of the NLL scheme 
and scale dependences and also discuss the sensitivity of our NLL-BFKL descriptions with respect to the photon 
impact factors. We also study the NLOQCD predictions, testing their relevance by comparing the use of different 
parton densities and different renormalization and factorization scales. 

The plan of the paper is the following. In section II, we present the phenomenological NLL-BFKL formulation 
of the forward-jet cross section for the two schemes called S3 and S4, while briefly highlighting the principles of its 
derivation. In section III, we compare the predictions of the two NLL-BFKL schemes with the data, and also with 
LL-BFKL and NLOQCD predictions. We discuss the dependence of our results on the choice of the hard scale with 
which a s is running in Section IV, and on the assumption made for the photon impact factors in Section V. Section VI 
presents the scale and parton-density dependences of the NLOQCD predictions. Section VII is devoted to conclusions 
and an outlook. 



II. FORWARD-JET PRODUCTION IN THE BFKL FRAMEWORK 



Forward-jet production in a lepton-proton collision is represented in Fig.l with the different kinematic variables. 
We denote by t/s the total energy of the lepton-proton collision and by Q 2 the virtuality of the intermediate photon 
that undergoes the hadronic interaction. We shall use the usual kinematic variables of deep inelastic scattering: 
x = Q 2 / (Q 2 + W 2 ) and y = Q 2 /(xs) where W is the center-of-mass energy of the photon-proton collision. In addition, 
kx 3> A.qcd is the transverse momentum of the jet and xj its longitudinal momentum fraction with respect to the 
proton. The QCD cross section for forward-jet production reads 

y + ^r) it ,».2 + (i - y) , .,.2 ) > (!) 



dxdQ 2 dx j dk 2 -, nxQ 2 \\ 2 J dxjdk^ dxjdk\ 



where da^*f^ JX jdxjdk\ is the cross section for forward-jet production in the collision of the transversely (T) or 
longitudinally (L) polarized virtual photon with the target proton. 

In the following, we consider the high-energy regime x<^ \ in which the rapidity interval Y = log(xj/x) is assumed 
to be very large. The NLL-BFKL forward-jet cross section of our analysis is given by: 



da^p rX = as(k 2 T )a s (Q 2 ) t , 2 , f d 1 ( Q 2 ^ 



feffixj, k 2 T ) J ±L ^ # iL ( 7 ) e ^xQ)xe„b,^Q)]r (2) 



with the complex integral running along the imaginary axis from 1/2— ioo to 1/2+ioo. The running coupling is given 
by 

a(k 2 )=a s (k 2 )N c /ir=[blog(k 2 /A 2 QCD )]- 1 , with b = ^J^ 1 ■ (3) 



In formula {2}, the NLL-BFKL effects are phenomenologically taken into account by the effective kernel Xeff(j,ct). 
Let us now give further details on this approximation. 




FIG. 1: Production of a forward jet in a lepton-proton collision. The kinematic variables of the problem are displayed. Q 2 is 
the virtuality of the photon that undergoes the hadronic interaction, kr is the transverse momentum of the forward jet and 
xj is its longitudinal momentum fraction with respect to the incident proton. Y is the (large) rapidity interval between the 
two hard probes. 



The scheme-dependent NLL-BFKL kernels provided by the regularization procedure xnll{i,^>) depend on 7, the 
Mellin variable conjugate to Q 2 fk\ and u>, the Mellin variable conjugate to W 2 /Qkr- In this work we shall consider 
the S3 and S4 schemes [H| . recalled in Appendix A, in which xnll is supplemented by an explicit a dependence. 
One writes the following consistency condition [ill [l7| 

w = a xnll (4) 

which represents the diagonalized form of the NLL-BFKL evolution equation and allows one to formulate the cross 
section ^ in terms of x e //(7,a). The approximation amounts to introduce the effective kernel to satisfy the con- 
sistency condition. Indeed, the effective kernel Xef f{l is defined from the NLL kernel Xnll{i,u) by solving the 
implicit equation 

Xef f (7, a) = xnll [if, a Xef f (7, a)] , (5) 

as a solution of the consistency condition ((3]). 

To highlight how the effective kernel enters in the formulation of the forward-jet cross section, let us consider the 
following inverse Mellin transformation over oj, the variable conjugate to e Y (~ the energy squared), where 7(7,0;) 
represents next-to-leading order corrections to the LO impact factors cp : 

duj I(7,u;)e" Y = 7[ 7 , axeffbf, a)] e « W/ ( 7 ,a)y ^ with ^ = d X NLL (g) 

2inuj - axNLL(l,u) 1 - ±nll[i, aXeffil, a)] du 

The factor in front of the exponential is an unknown correction, due both to the yet unknown next-to-leading order 
corrections to the LO impact factors and to the approximations made in satisfying the consistency equation through 
the effective kernel method. For simplicity, in the 7 integration of P]). we choose to factor this term out and treat it 
as a constant normalization parameter. 
Some other comments are in order. 

• In formula the renormalization scale is k 2 — kxQ, in agreement with the energy scale (l8j . In practice, 
one solves ([S]) with a — a(fcxQ)- Therefore, to each renormalization scale corresponds an effective kernel [l5j . 
In Section IV, we shall test the sensitivity of our results when using k 2 — A kxQ and varying A. Following 
formula ([5|), the effective kernel is modified accordingly for each scheme, and we also modify the energy scale 
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FIG. 2: Comparison between the NLL-BFKL effective kernels Xe//(7j«) obtained by solving the implicit equation © in the 
S3 and S4 schemes. The left plot shows Xe// (7, 0.25) as a function of 7 while the right plot shows the minimum Xe//(0.5,a) 
as a function of a. The dotted curves show the LL-BFKL kernel xll{j)- 



As we pointed out already, in formula (J2j) we use the leading-order (Mellin-transformed) impact factors 

4>Mi)\-„„ M 2^ 2 1 r 3 (i + 7 )r 3 (i- 7 ) ^(i + 7) (2- 7 ) 



4>l(l) J ~™ em ^ c ^ e «2 7 2 r(2-2 7 )r(2 + 2 7 )(3-2 7 )(, 2 7 (1 - 7) J' (?) 

for a transversely (T) and longitudinally (L) polarized virtual photon where et is the charge of the quark 
with flavor /. We consider massless quarks and sum over four flavors in l[7|). This is justified considering the 
rather high values of the photon virtuality (Q 2 > 5 GeV 2 ) used for the measurement. We point out that our 
phcnomcnological approach can be adapted to full NLL accuracy, once the next-to-leading order impact factors 
are available (the jet impact factors are known at next-to-leading order [19]). For completeness, we shall discuss 
the sensitivity of our results to typical next-leading modifications of (ftp ^(7) in Section V. 

In formula (0), f e ff{xj,kp) is the effective parton distribution function and resums the leading logarithms 
log(fc|7 Aq CD ). It obeys the following expression 

C 

feff(xj,kp)=g(xj,kp) + —(q(x J ,k' 2 ,) + q(xj,kp)) , (8) 

where g (resp. q, q) is the gluon (resp. quark, antiquark) distribution function in the incident proton. Since the 
forward-jet measurement involves perturbative values of kx and moderate values of x,j, formula (J2J) features the 
collinear factorization of feff, with k% chosen as the factorization scale. 

By comparison, the LL-BFKL formula is formally the same as with the substitutions 

Xeff — ► Xll{i) = 2-0(1) - "0(1 - 7) - "0(7) 7 a{k 2 ) — > a — const, parameter , (9) 

where r tp("-/) = d\ogT( , y)/d'y is the logarithmic derivative of the Gamma function. One obtains 



dxjdk 2 ~ k%Q 2 WW'KtJJ 2in\k? 



feff(*J,k 2 T ) £ J f T Al) e aXLLhW (10) 



Inserting formula ^ (resp. (|10p) into ([T]) gives the forward-jet cross section in the NLL-BFKL (resp. LL-BFKL) 
energy regime. In the LL-BFKL case, this is a 2-parameter formula: the overall normalization and a. In the NLL- 
BFKL case, each set of scales (Q 2 , k^) defines the running coupling constant and therefore we deal with only one free 
parameter, the overall normalization. The interesting property of our phenomenological approach is that formula (J2|) 
has formally the structure of the LL formula, but with only one free parameter and a NLL kernel. The delicate aspect 
of the problem comes from the scheme-dependent effective kernel Xef f ■ A comparison between the LL-BFKL kernel 
and the S3 and S4 effective kernels is shown in Fig. 2. As is well known, the NLL modifications to the BFKL kernel 
are quite important and will play an important phenomenological role in our analysis. 



III. NLL DESCRIPTION OF THE HI DATA 



The NLL-BFKL formula for the fully differential forward-jet cross section is obtained from (fTJ) and @. To compare 
the corresponding prediction with the data, one has to carry out a number of integrations over the kinematic variables. 
They have to be done while properly taking into account the kinematic cuts applied for the different measurements. 
The procedure is the same as the one described in Ref.[9j|, Appendix A. First one chooses the variables that lead 
to the weakest possible dependence of the differential cross section (we noticed that the best choice is 1/Q 2 , 
log(l/a;j), and log(l/x)) and then the integrations are computed numerically following the experimental cuts defined 
in 0.0;. 

To fix the normalization (the only free parameter) and check the quality of the data description using the BFKL 
formalism, we start by fitting the da/dx HI data The choice of this data set corresponds to the kinematical domain 
where the BFKL formalism is expected to hold (i<1 and Q 2 /k\ ~ 1). We then use the relative normalizations 
obtained between the different NLL BFKL calculations (S3 and S4) to make predictions for the triple differential 
cross section da / dxdk\dQ 2 . For this first analysis, the coupling a is running with the scale kxQ- 

A. The cross section da/dx 

We considered two kinds of fits: the first one is performed using statistical and systematics errors and the second 
one with statistical errors only. The systematics errors are very much point-to-point correlated, and this is why it 
is important to perform the fits with statistical errors only. Ideally, one should use the statistical errors added in 
quadrature with the uncorrelated ones but this information is not available. 

The fit results to the da/dx HI data are given in Table I. The \ 2 values (per degree of freedom) of the fits performed 
using the full (statistical and systematics) errors are quite good, always less than 1, for the two NLL-BFKL schemes we 
considered. This shows the possibility of describing the forward-jet cross section using the BFKL formalism at next- 
to-leading logarithmic accuracy. The fits using statistical errors only (which assume implicitly that the systematics 
are maximally correlated which is close to reality) bring about more constraints and show interesting features. The S4 
fit can describe the data better (x 2 = 10.0/5 d.o.f.) whereas the S3 scheme shows a higher value of % 2 (x 2 = 29.5/5). 
This indicates that the S4 scheme is favored. 

The curves corresponding to the fit with statistical errors only are displayed in Fig. 3 (upper plot), and they are 
compared with the LL-BFKL results taken from [§]. We notice the tiny difference between the LL and NLL results 
(the corresponding curves are barely distinguishable on the figure). This confirms that the data are consistent with 
the BFKL enhancement towards small values of x. Contrary to the proton structure function F2, the forward-jet cross 
section da/dx does not show large NLL-BFKL corrections, once the overall normalization fitted. This is due to the 
rather small value of the coupling a ~ 0.16 obtained in the LL-BFKL fit [9j, corresponding to an unphysically large 
effective scale. 

We also present in Fig. 3 the fixed order QCD calculation based on the DGLAP evolution of parton densities. The 
next-to-leading order (NLO) prediction of forward-jet cross sections is obtained using the NLOJET++ generator [20| . 
CTEQ6.1M [2l| parton densities were used, and the renormahzation [i r and the factorization scale [if were set equal 
to n 2 — [i 2 = Qk™ ax , where k™ ax corresponds to the maximal transverse momentum of forward jets in the event. 
The NLOQCD predictions do not describe the data at small values of x, as they are lower by a factor of order 2. 
The sensitivity of these predictions to variations of the renormahzation and factorization scales will be discussed in 
Section VI, as well as their dependence when using other parton densities. 

The fit parameters obtained with statistical error only will be used in the following to make predictions for other 
observables, namely the triple differential cross section da / dxdk^dQ 2 . The value of a for the LL-BFKL fit will be 
kept as well as the normalizations of the different NLL-BFKL calculations (S3 and S4). 



scheme 


fit 


X 2 /dof 


N 


S3 


stat. + syst. 


1.15/5 


1.591 ± 0.089 


S3 


stat. only 


29.5/5 


1.640 ± 0.019 ± 0.204 


S4 


stat. + syst. 


0.48/5 


1.356 ± 0.076 


S4 


stat. only 


10.0/5 


1.374 ± 0.016 ± 0.172 



TABLE I: Results of the NLL-BFKL fits to the HI da/dx data. The relative values between the different overall normalizations 
N can be compared. 
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FIG. 3: The forward-jet cross sections da/dx (upper plot) and da j dxdQ 2 dk% (lower plot, in nb/GeV 4 ) measured by the HI 
collaboration. Comparison with the two NLL-BFKL parametrizations S4 and S3 using the krQ scale, and with the LL-BFKL 
and NLOQCD predictions. In the case of da/dx, we see a good agreement between the data and the BFKL fits (the NLL-BFKL 
fits and the LL-BFKL fit are barely distinguishable on the figure) while the NLOQCD predictions do not describe the data. 
For da / 'dxdQ 2 'dk% , the best description of the data over the whole kinematic range is obtained in the NLL-BFKL approach. 
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FIG. 4: Renormalization scale dependence for the da/dx cross section. Upper plots: S3 scheme, lower plots: S4 scheme. The 
left plots show the results of the 2krQ and krQ/2 scales. The right plots show that in the considered range of renormalization 
scale, deviations from the scale krQ essentially affect the overall normalization. Therefore the quality of the fits (which only 
adjust the normalization) discussed in Section III is not altered by renormalization scale uncertainties. 



B. The triple differential cross section da/dxdk^dQ 2 

The triple differential cross section do / dxdk^dQ 2 is an interesting observable as it has been measured with three 
different fc|, and Q 2 cuts, yielding nine different regions for the ratio r — k\jQ 2 . It was noticed in that the LL-BFKL 
formalism leads to a good description of the data when r is close to 1 and deviates from the data when r is away from 
1, as effects due to the ordering between Q and fcy start to set in. NLOQCD predictions show the reverse trend. 

The HI data for do j dxdk\dQ 2 are shown in Fig. 3 (lower plot), and they are compared with the S4 and S3 
predictions, the LL-BFKL results (taken directly from 9]) and NLOQCD calculations. It is quite remarkable that the 
NLL-BFKL calculation, which includes some ordering between Q and kx, leads to a good description of the HI data 
over the full range. As was the case for da/dx, the difference between the LL and NLL results is small when r~ 1. 
By contrast when r differs from 1, the difference is significant, and the observable do / dxdk^dQ 2 is quite sensitive to 
NLL-BFKL effects. As a result, the best overall description of the data is obtained with the NLL-BFKL formalism. 



IV. RENORMALIZATION SCALE DEPENDENCE OF THE NLL DESCRIPTION 



In this section, we study the renormalization-scale dependence of the NLL-BFKL description. In the previous 
section, the choice was k 2 = kxQ and we now test the sensitivity of our results when u sing k 2 — A kxQ with A = 2, 
A = l/2, \ = kx/Q, or A=(Q//ct). We use formula ^) with the appropriate substitution [22| 

a{k T Q)^a(Xk T Q) + b a 2 (k T Q) log(A) (11) 

and with the effective kernel modified accordingly following formula ©• We also modify the energy scale kxQ^XkxQ- 
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FIG. 5: Renormalization-scale dependence of the triple differential cross section da / dxdk^dQ 2 (in nb/GeV 4 ) for the S4 scheme. 



We first study the case of da/dx and the results for the S3 and S4 schemes are shown in Fig. 4. We also display 
the results in terms of ratios with the prediction of the krQ scale chosen as the reference. We notice that the change 
of scale essentially affects the overall normalization and thus does not alter the quality of the fit, after readjusting 
the normalization. This is confirmed by the results of Table II (left table), which presents new fits performed to the 
da/dx data for the two scales 2 krQ and krQ/2 : for each scheme, the % 2 values are almost insensitive to the scale. 
Note that, for larger values of A such as A = 4, the quality of the fit deteriorates due to the large values reached by 
q in the non-perturbative range. Finally, due to the cut 0.5 < k\/Q 2 < 5 used for the measurement @, using the 
scales fcf. or Q 2 (instead of kxQ) yields much less modifications and changes in x 2 - 

The next step is to study the effect of the renormalization-scale dependence on the triple differential cross section. 
Applying the normalizations of Table II while computing the corresponding predictions gives, for either scale, results 
quite similar to those of Fig. 3. In Fig. 5, we rather show the predictions of the S4 scheme with the common normaliza- 
tion of Table I, and the scale dependence is generally found to be small. The biggest effects are uncertainties of about 
the same magnitude as the experimental errors, and they are obtained for large values of r — k\/Q 2 . The conclusion 
is identical in the case of the S3 scheme. 



Scheme 


Scale 


X 2 /dof 


N 




Scheme 


Impact factor 


X 2 /dof 


N 




k T Q 


10.0/5 


1.374 ± 0.016 








10.0/5 


1.374 ± 0.016 


S4 


2k T Q 


9.8/5 


1.644 ± 0.019 




S4 


7 (l/2) 


82.6/5 


0.655 ± 0.008 




k T Q/2 


8.8/5 


1.118 ± 0.013 








23.0/5 


2.694 ± 0.038 



TABLE II: Impact on the fits to the HI da/dx data for the S4 scheme when using different renormalization scales (left table) 
or modified impact factors (right table). 




FIG. 6: Impact factor dependence for the da/dx cross section in the S4 scheme. The left plot shows the results of the 
fits. The right plot shows that in the exact-gluon-kinematics case, the change of impact factors essentially affects the overall 
normalization, except at very small x, and thus modifies only slightly the quality of the fit. By contrast in the 7 (l/2) case, 
the \ 2 values are sensitive to modifications of the impact factor. 



V. IMPACT FACTOR DEPENDENCE OF THE NLL DESCRIPTION 



In this section, we study the dependence of our results on typical next-leading modifications of the impact factors. 
Indeed, since the next-leading impact factors are unknown, it is useful to see if our results are stable for different 
possible modifications. From |T]) and 1(2} one notes that the impact factors are involved in the computation of the 
fully differential cross section via the following factor in the integrand: 



(I- 2 / + y 2 /2)^(7) + (l-2/)0i(7) 



(12) 



In the previous section, the impact factors (ftp and <jP L were computed as functions of 7, and we treated unknown NLO 
corrections as a constant parameter. We now test the sensitivity of our results when using other prescriptions. 

Our first choice is to compute the LO impact factors at 7 = 1/2 and factor them out of the 7 integration as well. 
Another prescription that we shall study is the following: implement the higher-order corrections into the impact 
factor that are due to the exact gluon kinematics in the 7* — ► qqg transition. These have been calculated in [23| and 
can be taken into account in our approach: this is done using the lo— dependent Mellin-transformed impact factors 
02 gfe (7,w), where we recall that ui is the Mellin variable conjugate to W 2 /Qkr- Denoting (5 = w — 27+I, one writes: 



0^(7, oj) 



ira e 



r(7+<5)r( 7 ) 



' 9 47 2 r(^)(4-<5 2 )(5 2 -l) 



3(w + l) 2 + <5 2 + 8 
24(7+5)7 



^(7+^-^(7) / lu 2 [3(uj + 1) 2 + 9] + {5 2 -1)(S 2 - 
2lo8 \ 8(7 + 5)7(3cj 2 -<5 2 + 1) 



2w-9) 



(13) 



In practice, following ([6]), these impact factors are evaluated at uj = a{kTQ)Xeff[l, aferQ)] prior to the 7 integration. 
This prescription is motivated by the fact that, in the proton structure function analysis, these higher-order corrections 
allow for an improved DGLAP analysis [24| : indeed, they match the fixed-order results (at NLO and NNLO) for the 
splitting and coefficient functions of the DGLAP approach. 

We first study the case of da/dx. Table II (right table) presents new fits performed to the da/dx data for the 
different prescriptions. The fit results are shown in Fig. 6 (left plot) and we also display them in terms of ratios 
with the prediction of the (p{j) case chosen as the reference (right plot). In the 0^(7, w ) case : we notice that the 
change of impact factors essentially affects the overall normalization, except at very small values of x where the gluon 
kinematics become more restrictive. As a result the quality of the fit, which readjusts the normalization, is only 
slightly modified. By contrast in the </> 7 (l/2) case, the shape is clearly different, which yields a bad \ 2 - 

After applying the normalizations of Table II (right table), we now compare the effect of the different impact 
factors on the triple differential cross section. The results are given in Fig. 7. The effect is found to be small in the 
exact gluon- kinematics case, and even slightly improves the description. By contrast, differences are important in the 
</> 7 (l/2) case, especially at large values of r = k 2 -,/Q 2 . 
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FIG. 7: Impact factor dependence of the triple differential cross section da/dxdk^dQ 2 (in nb/GeV 4 ) in the S4 scheme. The 
relative normalizations coming from the da/dx fits have been applied. 



VI. PDF AND SCALE DEPENDENCE OF THE NLOQCD PREDICTIONS 



To obtain the NLOQCD predictions of the forward-jet data, the photon-parton hard cross section is computed at 
next-to- leading order in a s , and the leading and next-leading logarithms a™ log" Q 2 and a™ log™ -1 Q 2 are resummed 
using the DGLAP equations Q which govern the evolution of the parton distribution functions (PDFs) of the proton. 

The prediction for the forward-jet cross section at next-to-leading order is calculated using NLOJET++ [2fj]. The 
renormalization scale n 2 and the factorization scale /j, 2 were chosen as /a 2 . = [ij = Qk™ ax , where k™ ax denotes the 
maximal transverse momentum of forward jets in the event. To obtain the uncertainty associated with the choice 
of the scale, we varied the scales in the conventional range Qk™ ax /A < /i 2 — /i 2 < AQk™ ax . Another scale choice 
fi 2 = [i 2 — k™ ax2 was also tested; however, it yielded a result located within the mentioned scale uncertainty bounds 
and thus shall not be considered further. We point out that the scale uncertainty is rather large in the low— a; regime. 
This is due to the large NLO correction suggesting that higher-order corrections might be significant. 

An additional contribution to the total uncertainty of the calculation is due to the choice of the PDFs of the proton. 
Throughout the calculation, we use the CTEQ6.1M PDF parametrization which provides not only the central PDF 
So corresponding to the best fit to data, but also 40 additional distribution functions , , i = 1 ... 20 devoted to 
uncertainty studies 2l|. The total PDF uncertainty AX of the observable X is then computed as 
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(AX) 2 = £ 



x(s+) - x(sr) 



(14) 



We noticed that the main contribution to the PDF uncertainty comes from the gluon PDF. 

The DGLAP calculation of the da/dx distribution measured by HI is presented in Fig. 8, upper plot. This approach 
clearly fails to describe the data for the low values of x. In comparison with Fig. 3, we see that da/dx is more sensitive 
to the BFKL dynamics than the triple differential distribution as the deviation from the data is obvious. 
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FIG. 8: Upper plot: scale and PDF uncertainties of the da/dx distribution; the particular contributions of 40 different 
PDFs from CTEQ6.1M to the total PDF uncertainty (yellow area) is depicted in the inset. Lower plot: scale and PDF 
uncertainties for the triple differential distribution do jdxdQ 2 dk\ (in nb/GeV 4 ). In both cases, the scales varying in the range 



0.25 Qk™ ax < (j% = $ < 4 Qk™ ax axe displayed (shaded gray area) 



The PDF uncertainty study is shown in the inset. The solid lines represent the ratio of the cross section calculated 
with the various PDFs , S~ , normalized to the cross section calculated with the central PDF Sq. The gluon PDF 
(i=15) has the greatest impact on the uncertainty and the other PDFs can be neglected. The uncertainty study of 
da I dxdk^dQ 2 is shown in Fig. 8, lower plot. The PDF uncertainty was obtained taking into account the gluon density 
(which practically dictates the overall PDF uncertainty) only. 

The main conclusion that can be drawn from Fig. 8 is that at low values of x, the NLOQCD results suffer from large 
uncertainties, and they indicate that NNLO calculations arc needed to obtain genuine predictions in this framework. 

VII. CONCLUSIONS 

We performed a phenomenological analysis of the HI forward-jet data, looking for effects of next-to-leading loga- 
rithms in the BFKL approach. Let us briefly summarize our results. 

• For the cross section da/dx, measured in the kinematical regime Q 2 /k\ ~ 1, we obtain a good description of 
the HI data by NLL-BFKL predictions (see Table I and Fig. 3, upper plot). In addition, the difference between 
the LL-BFKL and NLL-BFKL descriptions is very small once the overall normalization is fit. This confirms the 
validity of the BFKL description of [1, Q previously obtained with the LL formula and a rather small effective 
coupling. 

• In the case of the triple differential cross section da / dxdk^dQ 2 , the same conclusions holds when r = k 2 -,/Q 2 

In addition when r differs from 1, the NLL-BFKL description is quite different from the LL-BFKL one, as it 
is closer to the NLOQCD calculation. As a result, the best overall description of the data for da / dxdk\dQ 2 is 
obtained with the NLL-BFKL formalism (see Fig. 3, lower plot). 

• The renormalization-scale dependence of our results has been thoroughly studied and we showed the stability 
of the NLL-BFKL approach when using the scales kxQ, 2kxQ and For da/dx, the change of scale 
essentially affects the overall normalization and does not alter the quality of the fits (see Fig. 4 and Table II, 
left table). For da/dxdk\dQ 2 , the biggest effect yields an uncertainty of about the same magnitude as the 
experimental errors (see Fig. 5). 

• We want to stress that the HERA data allow for a detailed study of the NLL-BFKL approach and of the QCD 
dynamics of forward jets. In particular, it has the potential to address the question of the remaining ambiguity 
corresponding to the dependence on the specific regularization scheme of the NLL kernel. For instance, the 
predictions of the S3 scheme do not compare with the data as well as the predictions of the S4 scheme, as 
indicated by the x 2 values given in Table I. Therefore, it would be very interesting to compare the data with 
other regularization procedures [H|,[I2] than those used here. However, these other solutions proposed to remove 
the spurious singularities of the NLL kernel are not in such a suitable form for phenomenology, hence this issue 
will be addressed in a separate work. 

• Our analysis is to be completed with the next-leading photon impact factors, when available. However we tested 
the stability of our approach when implementing typical next-leading modifications of leading-order impact 
factors. The results show some sensitivity (see Fig. 6-7 and Table II, right table) with the 7 (l/2) prescription. 
It is however interesting that using the impact factors (|13|) with exact gluon kinematics (which, in the structure 
function case, allow for an improved DGLAP analysis) still gives a good fit of the da/dx data, and also a better 
description of the da / dxdk^dQ 2 cross section. This indicates that, when the next-leading impact factors will 
be known, our predictions could be stable. 

• Finally, we computed the NLOQCD predictions using the NLOJET++ generator and CTEQ6.1M [H| partem 
densities. We tested their relevance by comparing the use of different parton densities and renormalization and 
factorization scales (see Fig. 8). The NLOQCD predictions do not describe the data at small values of x, and 
they suffer from large uncertainties, showing the need for higher-order corrections in this framework. 

Forward-jet production is the first observable for which the NLL-BFKL description works, while the standard 
NLOQCD does not work. We need complete knowledge of the next-to-leading impact factors before drawing final 
conclusions, but our analysis strongly suggests that the data show the BFKL enhancement at small values of x. This 
is of great interest in view of the LHC, where similar QCD dynamics will be tested with Mueller-Navelet jets pM |26|. 
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APPENDIX I: The S3 and S4 regularization schemes 



In this appendix, we recall the regularization procedure of to obtain xnll{i, in the S3 and S4 schemes, 
starting point is the scale invariant (and 7 «-> 1 — 7 symmetric) part of the NLL-BFKL kernel 



The 



Xi(7)= 2^(3) + 



-0(7) 



7T 2 C0S(7T7) 



4shr(vT7)(l - 27) 



o 1 (, 1 N f\ 2 + 37(1-7) 
' : ' + 7V3;(3-2 7 )(l + 2 7 ) 



with b given in xll given in ((9]), and 

fc+1/2 



fc=0 



fc + 2 



fc + 1 



The pole structure of Xlil) at 7 = (and by symmetry at 7 = 1) is: 

di 



Xi{l) = - T - i +^ + - + 0{l) 



with 



N f ( 13 \ 

d\ = — 5 H - 

18N C \ 2N* J 



11 N f / 2 

d 2 = + — — 1 » 

8 12N r \ N? 



• The S3 scheme kernel XS3{j,u) is given by 
XS3(l,u) = [1~ aA 



i , . 2aB + Lu\ / 2aB + uj 
2V>(1) -V) 7+ 5 -V 1-7+ = 



+a ^ xi (7) + ^Xll(7) + ( B + ) [^(7) + V>'(1 - 7)] 



with A and -B chosen to cancel the singularities of Xi{l) at 7 = : A = —d\ — tt 2 /6 and B — —d%. 
• The S4 scheme kernel xsa{i,u) is defined with the help of the function 7(7) = 1/7+ 1/(1 — 7) : 



XS4(7,w) = XLL{l)-f{l) + [l-aA]f{u + 2aB, 1 ) + a\xi(H)+Af{ 1 )+ [B + 



In this scheme, A and B are given by A = — d\ — 1/2 and B = —d 2 - 



XLL (7) 



1 



(15) 



(16) 



(17) 



(18) 



(19) 



7 2 (1-7) 5 



(20) 



APPENDIX II: Comparison between the exact NLL-BFKL 7 integration and a saddle-point approximation 



It is possible to estimate the complex integration in @ using a saddle-point approximation in 7. In the BFKL 
regime we are working in, Y is very large, and the saddle-point equation 

log^/Q 2 ) 



dXeff 



(7c a) = Xe//(7c,a) 



aY 



(21) 



becomes x' e f f (7ci a ) — 0- Hence one finds for the theoretical forward-jet cross section 



dal*V JX a s (fc 2 H(Q 2 ) 



T.L 

dxjdk 2 -, 



/ I f ff ( Tj ( 



f T ,M /_ log 2 (Q 2 /^) 

exp ax e //(7c,a)r 



27rax"//(7c,a) V 



where x'/ ff ^d 2 Xeff/d-f 2 



2cc X 'l f f{lc,cc) Y J ' 
(22) 
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FIG. 9: The forward-jet cross sections da/dx (upper plot) and da/dxdQ 2 dk^ (lower plot, in nb/GeV 4 ) measured by the HI 
collaboration and compared with two NLL-BFKL predictions of the S4 scheme, obtained with the exact 7 integration ^ or 
the saddle-point approximation (|22|l . The LL-BFKL and NLOQCD results are also recalled. 



For each set of scales (Q 2 , kj), it is possible to extract the values of j c , otXeff(lc, o>) and ax"//(7o 01) after solving 
the implicit equation (]5|). 

This approach was considered in [lfjj, and compared with the results obtained with the exact integration. In 
this appendix we discuss this comparison in more detail. At the level of the differential cross sections @ and p!?]) 
there are some differences between the exact NLL-BFKL integration and the saddle-point approximation. But when 
considering the integrated, experimentally measured cross sections da/dx and da / dxdk^dQ 2 , the description of the 
data is similar. This is shown in Fig. 9. 

We recall that, starting from ([T]), one has to carry out a number of integrations over the kinematic variables, 
which have to be done while properly taking into account the kinematic cuts applied for the different measurements 
Q. This procedure seems to erase the differences between the exact NLL-BFKL integration and the saddle-point 
approximation. This is even more so for da/dx. 
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